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We investigate systems of real scalar fields in bidimensional spacetime, dealing with potentials 
that are small modifications of potentials that admit supersymmetric extensions. The modifications 
are controlled by a real parameter, which allows implementing a perturbation procedure when such 
parameter is small. The procedure allows obtaining the energy and topological charge in closed 
forms, up to first order in the parameter. We illustrate the procedure with some examples. In 
particular, we show how to remove the degeneracy in energy for the one-field and the two-field 
solutions that appear in a model of two real scalar fields. 



I. INTRODUCTION 



Domain walls are defect structures that appear in systems engendering spontaneous breaking of discrete symmetry. 
They are of interest for instance in Condensed Matter, as interfaces in magnetic materials [y, as seeds for pattern 
formation , and as interfaces in ferroelectric crystals , and in Cosmology, as seeds for the formation of structures 
1^,0 in the early universe. Domain walls spring through the immersion of kink-like solutions of (1, 1) space-time 
dimensions to higher spatial dimensions. Standardly, kinks or domain walls appear in scalar field models, but they 
may also be present in extended systems, that include fermionic fields, with or without supersymmetry. 

Recently, the study of supersymmetric models has brought new issues, as for instance in the investigations of a Wess- 
Zumino model |^,^ engendering the symmetry. This model allows the identification of a Bogomol'nyi equation |Tc| ] 
for a triple junction that breaks 1/4 supersymmetry of the model. We learn from these works that supersymmetry 
helps easing calculations concerning the presence of the triple junction 12 1. However, supersymmetry seems to 



play no central role when the issue is the tiling of the plane and not the triple junction ifself |13Hl5|. Furthermore, 
in Ref. it was shown how to entrap a planar regular hexagonal network of defects inside a domain wall. This 
model involves three real scalar fields, and engenders the Z2 x symmetry, but it seems to have no supersymmetric 
extension. 

These facts motivate the study of bosonic models that do not support supersymmetric extensions. However, to keep 
track of supersymmetry we investigate models that are close to the real bosonic portions of supersymmetric systems. 
These models are defined by potentials that contain two parts: the first part defines the real bosonic portion of a 
supersymmetric system; the second part defines the extended model. The first part alone constitutes the basic model, 
which supports static field configurations that minimize the energy, attained by static field configurations that solve 
first order differential equations. In particular, we examine models that support topological solutions that belong to 
the same topological sector, and are degenerate, having the very same energy. An example of this was investigated in 
in a model of two coupled scalar fields. This model has been recently extended to the case of several scalar fields 



in |lq ]. We explore the possibility of extending this system, in order to remove the degeneracy between the one-field 
and the two-field solutions. This investigation is of intrinsic interest, and may also help examining applications to 
Cosmology and to Condensed Matter. In Cosmology we recall the usual route and also the new possibility 

In Condensed Matter, we envisage other issues, concerning for instance the presence of Ising and Bloch wall 
interfaces in magnetic materials described by the anysotropic XY model |^,^ , and the structural phase transition in 
ferroelectric crystals |^,|l9|. 

The real bosonic sector of supersymmetric systems described by n chiral superfields $" contains n real 

scalar fields 0i, 02, The potential V is written in terms of the superpotential W, in a way such that for 
V = V{(pi,(p2, ■■■, 4>n) and W = 02, ■■■An) we get 

y = \wi + \wi + ■■■ + \wi (1) 

where W^. = dW/d(j)i,i = 1,2, ...,n. These systems have been investigated in several different contexts in ||l5| , [l7| , |l9| , 
and in references therein. In the present work we investigate systems where the potential includes an extra piece, that 
modifies the above potential according to 

V,=V+^eF (2) 
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where e is a parameter, real, and F = -F'((/>i, 4>2, (j^n) is in principle an arbitrary function of the fields. In the present 
work we show that if e is small, we can develop a perturbation procedure that gives closed results up to first order in 
e. The perturbation procedure is based on former investigations [|2^,^, and works nicely for potentials of the form 
when the function F obeys some restrictions. We start our investigations using natural units, but we work with 
dimensionless fields and coordinates, and sometimes we refer to the model with e = as the primary system, and to 
the complete model, with e 7^ 0, as the extended system. 



II. GENERAL CONSIDERATIONS 



In order to better understand the procedure, let us first consider models that are bosonic portions of supersymmetric 
theories. In this case we have 

n 

C = -Y,d»cb,d'^(b^-Vi(l)i,...,cj,n) (3) 

i=l 

where V is given by Eq. (|l|). We work in (1, 1) space-time dimensions. The equation of motion for static fields are 



dx'^ 



W^,W^,^^ + ... + M^0„W^0„02 



= W^,W^,^^ + ... + T4^0„W^0„0„ 
These second order equations are solved by field configurations that solve the first order differential equations 

^^M^, ^ = iy^ ^ = 
dx dx ■'■ dx 

These are the Bogomol'nyi equations. The energy of the static solutions can be written as 



E 

2 



1 pOC ^ 



dx 



2 

2 



(4) 



In general, the system may have several distinct sectors, which may be identified by the two vacuum states the static 
solution connect. Thus, if we use the set of numbers {wa,f(),Wc, •■•} to mark the vacuum states of the model, we can 
write 



E'*.Et.\\yti^l-w,:) (5) 



Here = | Wah|, and Wab = M^(wa) — W{vh). A pair of vacua defines a topological sector, and in the sector (a6) the 
energy is minimized to the bound E'^ for field configurations that obeys the above first order equations. Solutions of 
the Bogomol'nyi equations are named Bogomol'nyi-Prasad-Sommerfield solutions pO|j2^ , and we refer to this bound 
as the BPS bound. It is possible that W(vo) = W{v\j), giving a vanishing Wah- In this case the topological sector 
cannot support BPS states, and we refer to this as a non-BPS sector - see Ref. The topological features of these 
solutions can be accounted for by introducing for instance the topological current 



" = \ e'^^da'i (6) 



J 

where $ is a column vector, such that <I>* — (0i 02 ■ • • f^n)- 

We illustrate the general situation with some examples. We work with natural units. However, we avoid unimportant 
considerations by considering the Lagrangian density jC, written in terms of dimensionless fields and coordinates. The 
usual form of the Lagrangian density is CJ = 7/I, where 7 is a constant, which carries the correct dimension of the 
Lagrangian density in natural units. 
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A. One Real Scalar Field 



In the case of one real scalar field we have 



The potential V{(j)) speficies the system. As an example we consider the model 

= 102-101 + 1 (8) 

This potential was recently considered in Ref. There are two minima, at the values {vi — 1,V2 = — !}■ The 

equation of motion is (P(j)/dx'^ = <l> — 0/101- It has the solutions (j>±{x) = ±(a;/|a;|)[l — exp(— |a;|)]. This model can be 
written in terms of the function 

W^(0) = - 1 101 (9) 
This means that the topological sector is a BPS sector. The BPS solution satisfies 

^ = 1 - 101 (10) 
ax 

This solution is 0(x) = (x/|x|) [1 — exp(— |a:|)]. It has a kink-like profile, and is linearly stable minimizing the 
energy to Eb = 1. 

B. Two Real Scalar Fields 

We exemplify the case of two real scalar fields by considering the superpotential 

W^(0,x) = 0-l0'-r0x' (11) 
where r is a real parameter. Here the potencial gets to the form 

V{cb, x) ^\{^^- 1)' - -rx^ + Kl + 2r)02x' + \r'x^ (12) 

There are two minima for r < 0, wi = (1, 0), V2 = (—1, 0). For r > thare are four minima, the two former ones, and 
also the other two vg, = (0, y^l/r), = (0, —y^l/r). 
In this model the first order equations are 

(13) 

ax 

We see that for x ^ the first order equations demand dip/dx = (1 — 0^), and the BPS defect solution is 4>{x) = 
tanh(x). However, for 0-^0 the first order equations demand that — 1/r, that is, that the x field should be at 
the corresponding minima. These results are manifestations that the sector defined by the minima (±1,0) is a BPS 
sector, while the other sector, defined by the minima (0, ±^/T/r) is a non-BPS sector, and the non-BPS solutions are 



stable if and only if l/r < 1 |15 



The solution = tanh(2:) and x = is a one-field solution, but the system also supports two-field solutions. In the 
sector that connect the minima (±1,0), the two-field BPS solutions have the explicit form 

0(a;) = tanh(2ra;) (15) 

X(x) = ±\j- - 2 sech(2ra;) (16) 

which requires that l/r > 2. The one-field and two-field solutions that appear in the sector connecting the minima 
(±1,0) have the same energy, E = 4/3. In the (0, x) space the real line a; g R is mapped to line segments: for 
the one-field solution we get a straight line segment going from (—1,0) to (1,0), and for the two-field solutions we 
get an ellyptic arc, obeying 0^ + x^ /{l/r — 2) = 1. These solutions map the Ising and Bloch walls that appear in 
magnetic systems, respectively, and represent solutions of the anisotropic XY model, that describe interfaces between 
ferromagnetic domains P,p^. 
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III. EXTENDED SYSTEMS 



Here we deal with systems defined by extended potentials, that differ from the above primary potentials by some 
slight modifications. We start dealing with the general model 

1 " 1 

^^=2^^^'' + 2'^^'^^'-''^"^ ^^^^ 

i=l 

where e is a real parameter, infinitesimal. In this case we can use such parameter to control the procedure of extending 
the result beyond the e-independent value. The correction to the potential depends on F — F{4>i, ■■■,4'n)- In general, 
we may have two different types of functions: functions that respect and functions that do not respect the symmetry 
of the original system. Both type of functions are important to describe situations where the system is modified by 
the presence of external fields, chemical potentials, etc. However, in the present work we are interested mainly on 
topological solitons, in investigating the topological sectors of the model. Thus, we consider the case of functions 
that respect some symmetry of the primary system. In this case the function F accounts for modifications of the 
original system, without destroying the topological sector one is investigating. This means that in the set of possible 
vacuum states {va,VhTVc, ■■■} of the primary system, at least the vacua Va and Vh remain present, althought they 
can be slightly changed to v% and vl- We shall be investigating slight modifications in the topological sectors of 
the BPS type, modifications that do not destroy the sectors themselves. In this case we write the static solution 
4>\{x), 4>l2{x), 4>%{x) of the extended system in terms of the static solution of the original model in the form, up to 
first order in e 



<j)l{x) = (f)i{x) + erii{x), i = l,2, 
In the extended system, the energy of static solutions has the form 



2 



dx 



1=1 



\ dx 



wi 



1 



dx F{(j>i,(j>2, 



It can be written as, in the case of the BPS sector that connects the vacuum states labeled by a and b, 

2 



dx F{(j)i,4>2, ■■■,4'n) 



(18) 



(19) 



(20) 



We now follow Ref. ||2^ . We see that the field 0^ given by Eq. ( p^ ) shows that the first term in the above expression 
for the energy do not contribute to first order in e. This fact allows writing the energy in the form 



1 



dx F{4>i,4>2,...,(t>n)+0{e') 



(21) 



Here E1^ = \Wab\ and Wab = W{v'^^) — W{vl). The correction to the potential is small, and we separate two cases: 
first, the case where the correction does not change the minima of the primary potential, that is, the case where 
= Va and vl = Vb', and second, the case where the correction slightly modifies the minima of the primary potential. 
In the first case both (j^ti^) ^nd have the same asymptotic behavior, thus r]i{x) must vanish asymptotically; in 

the second case </>|(a;) is asymptotically different from 4>i{x), thus T]i{x) cannot vanish asymptotically. In both cases 
the term E'^ exactly reproduces the corresponding term in the primary system, up to first order in e. This is so 
because W[(l)\{±oo), ^^(±00)] can be expanded to give 



VK[^i(±c»),.. .,(/.„ (±00)]+, 



i=l 



dW 

rjt (±00) —— 



.(±oo) 



However, we know that (f)i{ztoo) are minima of the primary model, and are extrema of W. Thus, the second term 
in the above expression vanishes. For the topological charge, in the first case we see that it does not change, giving 
Qt = Qt- However, in the second case the topological charges change according to = Qt + e AQ, where Q^, Qt 
and AQ are n-component vectors, and AQ accounts for the difference r]i{oo) — r]i(— 00), i — 1,2, ...,n. 

The above results are general results, and we illustrate the general procedure with some examples, spliting the 
investigation in the two subsections that follows, which deal with one and two real scalar fields separately. 
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A. The case of one real scalar field 



In the case of one field, let us first consider the 0^ model, defined by the potential 

V{<t^) = \{<f'-lf (22) 

We consider Fi{(j)) = (0^ — I)^, which is an example where the correction does not change the minima of the primary 
system, and so = Qt = 1- In this case the energy becomes 

E = Eb (l + U) (23) 



where Eb = (4/3). It can be greater (e > 0) or lesser (e < 0) than the energy of the unperturbed system. 
We notice that the above Fi(0) allows rewriting the potential as 

n0)-^(i+e)(<A'-i)' (24) 

This potential requires that e > — 1, and shows that the extended system is very much like the primary model, with 
the coupling for self-interaction changed by the e term. In this case we can find the energy of the static solution 
exactly. It is 



E = EbVT+V (25) 



For e very small, we expand the above result to get the former answer, Eq. (23), and this shows that our approach 
works correctly. 

As a second example, let us consider the same primary model and another function 

^^2(0) = (1-0') (26) 

In this case the correction does change the minima of the primary potential. The new minima are at ±(1 + e/4). The 
energy is 

E = EB(^l + le^ (27) 

It can be greater (e > 0) or lesser (e < 0) than the energy of the unperturbed system. The topological charge changes 
to Q'^ = Qt{1 + €/4). 

Here we also notice that with the above correction of Eq. ( p6| ) the potential can be rewritten in the form 



which is correct to first order in e. In this case the energy of the static solution is 

1 \ 3/2 



E = EBi^l + -ej (29) 

However, since e very small, we expand the above result to get the former answer, Eq. (|2^). We notice that the 
modification in the minima of the potential changes the topological charge, and also the energy of the topological 
solution. The energy increases or decreases, depending on the increasing or decreasing of the spontaneous symmetry 
breaking parameter. 

We now consider the model defined in Eq. (||). We extend this model with the above F2 (</*)■ The new potential is 

vM-l^'-l'l^l + l + l^i-^') (30) 

For e 7^ 0, small, the minima change from v± — ±1 to v^. — ±(1 + e). The energy of the topological solution changes 
to 
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We see that the energy decreases when the spontaneous symmetry breaking parameter decreases. 
We notice that the above potential can be written as 



(31) 



This resuh is vahd up to first order in e. It aUows introducing a superpotential, and the first order equation is 

101 



It has the BPS solution 



ax vl + e 



0(x) = (l + e)^(l-e-l-l/(i+^)) 



(32) 



(33) 



(34) 



This shows that both the amplitude and width of the former kink-like solution change in the extended model. The 
energy is (1 + e)'^/^, but since e is small we can expand this result to get the former answer, given by Eq. (^l|). 



B. The case of two real scalar fields 



The case of two fields is more involved, and we envisage several distinct possibilities of illustrating this situation. 
We consider the example presented in Sec. ||. We extend that model with the function 



Fi(0,x) = (l-0') 

In this case, in the BPS sector defined by the minima (±1,0) the correction to the one-field solution is 



EI = Eb[i + -e 



In the case of the two-field solution (|T^ and (^6|) we get 



El= Eb\1 + -e 



3 1 



8 r 



We introduce the ratio between energies, R = E2/E1. We see that 



We consider another perturbation 



It gives 



and 



«.^l + 5.(l-2 



^2(0, x) = rx^ 



El — Eb 



E2 — Eb 



They give 



(35) 
(36) 

(37) 

(38) 

(39) 
(40) 

(41) 
(42) 
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We see that i?i = i?2- The ratio R = E2/E1 does not depend on the way one extends the model, using Fi = (1 — (fp') 
or F2 = ry^ . We notice that the extension with Fi = (1 — (fp') changes the minima of the primary model from (±1, 0) 
to (±1 ± e/4), 0). Thus, the topological charge of the sector changes according to 



1/4 




(43) 



The other extension, that uses F2 — rx^, does not change the minima (±1,0), so the topological charge in this BPS 
sector remains the same = Qt- The above examples show two distinct ways of removing the degeneracy between 
the standard one-field solution and two-field solutions (|l5|) and (|l6|). The two-field solutions are less energetic for 
e < 0. This means that one favors the non-trivial two-field configuration when the symmetry breaking parameter 
decreases in one or in the two </> and x directions. 

There is another way of removing the degeneracy of the one-field and the two-field solutions that appear in the 
sector connecting the minima (±1,0). We consider the case where the extra piece contains interactions between the 
two fields, for instance 



1,2, 



(44) 



In this case the minima of the primary system do not change, so the topological charge gets no modification. The 
same happens to the energy of the one-field solution. However, the energy of the two-field solutions changes to 



E'^^^ = Ef 



1 



2fc+ 1 



(45) 



We see that in the limit fc ^ ones obtains the former result, giving in Eq. (g^). We also see that the sign of e 
makes the energy density of the two-field solutions to be higher (e > 0) or lower (e < 0) than the energy density of 
the one-field solution, removing the degeneracy they have in the primary system. 

In the former model of two real scalar fields with the fimction Fi (</>,%) as in Eq. (^5|), the potential has the form 



(46) 



There are four minima, two at (0, x±), x± — ±-\/l/?', and two at {4>±, 0), 4>± = ±(l + e/4). We compare this potential 
with the potential of the primary model. The presence of the extra term shows that, for r > 1 the stable non-BPS 
solution implies 



1 



1) 



1 



red 



(47) 



and 



This shows that the (squared) mass of the </>— meson is, inside the kink 



m^(m) 



411 + 1. 



Outside the kink we get m^Aout) — Ar — e, which gives the ratio 



(48) 



(49) 



mliin) ^ 1 
mP^{out) r 



2r 



(50) 



Thus, if the non-BPS kink entraps 0— mesons in the primary system, the entrapment is still more efficient in the 
extended system, for e < 0. We see that deviations from the BPS bound may improve the efhciency of the mechanism 
for the entrapment of the other field. 
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C. Another case 



Let us now consider the model investigated in Ref. ||T^. It is described by three real scalar fields, and the potential 
has the form 

Via, </>, X) = ^ (a^ (ra' - ^) (0^ + x') + i^' + x'f - " ix\ (51) 

The projection with (i^, %) (0, 0) gives 

n-) = ^(^'-^)' (52) 

This potential can be written with the superpotential VF(i7) = (3-\/3/2)(T — (2V3/9)ct'^. This fact shows that the 
defect cr(x) = (3/2) tanh(\/3a:) that appears in this case is a BPS defect. We can then extend this model adding to 
the potential in Eq. ( ^l| ) a term depending on the a field, for instance 

We use former results to see that this term contributes to decreasing or increasing the energy of the basic defect, 
increasing or decreasing the efficiency of the mechanism for the entrapment of the network the other two fields (/) and 
X may generate. 



IV. COMMENTS 



In this paper we have examined systems described by real scalar fields, in which the energy of static field config- 
urations is in the vicinity of the BPS bound. This bound is attained by field configurations that solve first order 
equations, and minimize the energy. The BPS bound appears in the real bosonic sector of supersymmetric theories 
described by chiral superfields. The systems we have investigated are extensions of primary systems, described by 
potentials given by the functions W = Vl^((/'i, <i>2, ■■•) and F = -^(01, 02, ■••), the specific form 

+ \K. + • • • + + \ ^ Fiq>i,h, -An) 

The parameter e is real, and is used to control the deviation from the primary system, described in terms of the 
superpotential W. 

In the extended system, if the static solutions are similar to the solutions one finds in the primary system, that 
is, for (plix) = 4>iix) + erjiix), we could examine the energy related to the new defect solutions 4>lix), and write the 
first corrections in e in a closed form, independently of the specific form of the new defect solution itself. The formal 
results are of direct interest to field theory, where they may be used to improve the mechanism for the entrapment of 
the other field 

We have investigated specific systems, and we have found diverse possibilities of removing the degeneracy between 
different types of solutions, without destroying the degeneracy of the vacuum states. The examples we have pre- 
sented serve to illustrate some practical possibilities of removing defect degeneration, and this is of direct interest in 
application in specific physical situations. In ferroelectric crystals, for instance, the order parameters that control 
structural phase transitions may be changed by applications of external pressure along specific planar directions in 
the crystal. This is a tipical scenario for changing the parameters that control spontaneous symmetry breaking inside 
the crystalline material, changing the energetics of the structural phase transition. In the systems we have examined 
in the present work, the presence of external pressure may be directly mapped into specific forms of F, that controls 
the extended system. The present work opens a new route for exploring systems of coupled scalar fields, intending to 
mimic specific systems in applications to Cosmology and to Condensed Matter. We postpone to the near future some 
specific investigations. 
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